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Acoustic wave propagation in a fluid with a random assortment of identical cylindrical scatterers is 
considered. While the leading order correction to the effective wavenumber of the coherent wave is 
well established at dilute areal density (no) of scatterers, in this paper the higher order dependence 
of the coherent wavenumber on no is developed in several directions. Starting from the quasi- 
crystalline approximation (QCA) a consistent method is described for continuing the Linton and 
Martin formula, which is second order in no, to higher orders. Explicit formulas are provided for 
corrections to the effective wavenumber up to O(no). Then, using the QCA theory as a basis, 
generalized self consistent schemes are developed and compared with self consistent schemes using 
other dynamic effective medium theories. It is shown that the Linton and Martin formula provides 
a closed self-consistent scheme, unlike other approaches. 



PACS numbers: 43.35Bf, 43.20.Fn, 43.20.Hq 



I. INTRODUCTION 

It is often assumed that multiple scattering by a 
dilute array of scatterers in a perfect fluid may be 
described by the propagation of a coherent wave 
representing the acoustic field averaged over all 
possible scattcrer configurations. The coherent 
wave has a complex- valued wavenumber, fc e //, of- 
ten called the effective wavenumber. Its imaginary 
part accounts for loss due to scattering in all di- 
rections. Foldy first derived an expression for the 
effective wavenumber of point scatterers (isotropic 
scattering)^. Subsequently, Waterman and Tru- 
ell considered finite size scatterers (non-isotropic 
scattering)^ and obtained a second order correc- 
tion to Foldy's formula in terms of the scatterer 
density. The averaged exciting field in£ is obtained 
by replacing the spherical scatterers by point scat- 
terers with the same angle-dependent far-field scat- 
tering amplitude as the actual finite size scatterers. 
For cylindrical scatterers, Waterman and Trucll's 
approach provides k e f f ~ kwr whero^ 



k 2 



[*-^/(0)] 2 -[^/M] 2 , (i) 



with /(#), defined in ^JTlJ the far-field scattered am- 
plitude in direction 9 for each cylinder and no the 
number of scatterers per unit area. 

In 1967 Lloyd and Berryi proposed an explicit 
expression for spherical scatterers different from 
that of Waterman and Truell. While Lloyd and 
Berry's formula had been obtained for spherical 
scatterers by using a different method^, Linton and 
Martin derived its counterpart for cylindrical scat- 
terers in 2005^. Linton and Martin's formula may 
be recovered from Fikioris and Waterman's disper- 
sion relation^ by expanding it in powers of no, un- 
der the assumption that no/k 2 << 1, and by let- 
ting the radius of exclusion in the hole correction 
tends to zero. For cylindrical scatterers, Linton 
and Martin's approach yields fc e // = ^lai, 
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[f{9)] 2 A9. 
(2) 

This formula is only valid for symmetric scatter- 
ing functions satisfying f(9) = /(— 9)^. Its gen- 
eralization to arbitrary / is obtained by replac- 
ing f 2 (9) with f{9)f(-9) in ©, see below and£. 
Linton and Martin's formula is a consequence of 
the widely used closure assumption known as the 
quasi-crystalline approximation (QCA). Martin 
and Maurel^ demonstrated that QCA agrees with 
the Lippmann-Schwinger approach (weak scatter- 



Multiple scattering by cylinders 1 



ing) to 0(n§). 

The Independent Scattering Approximation 
(ISA)^ is a simplistic approximation in multiple- 
scattering theory, sometimes used without justi- 
fication. It can be deduced directly from Eq. (J} 
and Eq. ([2]) by neglecting the term of second order 
in n /k 2 , giving k e ff = k IS A with 

k 2 isA = k 2 - 4m /(0). (3) 

All of these methods arc explicit since k e ff is 
given by a formula. But the question remains: 
which method is the most accurate? This de- 
pends largely on the value of the concentration, 
c = no7ra 2 , where a is the radius for circular cylin- 
ders or the radius of the enclosing cylinder for non- 
circular scatterers. It is evident that ISA, which is 
a method of order one in no/k 2 = c/ir(ka) 2 , is the 
least precise. As for the two others, according to 
the analysis made by Linton and Martin in~, it is 
justified to think that k^M is more accurate than 
kwT- However, the difference between kwr and 
k^M is small at low concentration, as shown inii 
when c = 6%. Ref.— also shows that the Water- 
man and Truell method or ISA fail while Linton 
and Martin's gives better results for higher con- 
centration: c = 14%. 

There are two other methods that deal with the 
far field scattered amplitudes, but they are im- 
plicit; k e ff is obtained by solving an equation. The 
most famous is the Coherent Potential Approxima- 
tion (CPA), and its generalization (GCPA), both 
based on Dyson's equation^. There is also the 
Generalized Self Consistent Model (GSCM) which 
is derived using a self consistent scheme applied to 
Waterman and Truell's formula^. The existence 
of a non trivial solution to the homogeneous sys- 
tem of linear equations obtained by Fikioris and 
Waterman also allows one to calculate a k e ff^. 
According to the finding of Cowan et alJ^ the CPA 
is capable of accounting for concentrations up to 
60% for glass beads in a liquid mixture of water 
and glycerol. Of course, scatterers do not radiate 
as much at this level of concentration as compared 
to those considered at the lower levels of concen- 
tration ir*ii. Nevertheless, the results oJi£ suggest 
that implicit methods can be more powerful than 
explicit ones even if, to our knowledge, no rigorous 
comparison has been made between the two types 
of methods. Ref provides comprehensive numer- 
ical comparisons of eight different explicit and im- 
plicit methods. 

In this paper we search for explicit and implicit 
expressions of fc e // that could allow us to con- 
sider higher concentrations for multiple scatter- 



ing problems. For explicit methods, the only way 
to improve fc e // is to extend Linton and Martin's 
formula to orders higher than 2 in concentration. 
This is accomplished in iflll In this regard we note 
that Waterman and Truell's formula is of order 
2 in concentration, and that it is is not possible 
to go further. For implicit methods we use the 
self consistent scheme as used by Yang and Mai 
ir>l£. When applied to the ISA this self consis- 
tent scheme leads to CPA. It is applied to the gen- 
eralized Linton and Martin formula in £|IVI where 
the physical meaning of the self consistent scheme 
is discussed. In particular, we obtain a new re- 
sult that generalizes the CPA. In SjV] we compare 
all the effective wavenumbers in the Rayleigh limit 
(low frequencies). 

II. HIGHER ORDER THEORY 
A. Multiple scattering formulation 

The problem is formulated in terms of the pres- 
sure ijj(r) which satisfies the Helmholtz equation 
in the interstitial space between scatterers: 

V 2 V + k 2 ip = 0, (4) 

where k = lu/c and c is the speed of sound. Time 
harmonic dependence e~ tult is assumed. Consider 
first a system of N scatterers with fixed positions, 
for which the total field can be expressed 

N 

W) = Vw(r) + T(r^ E (r, ?,). (5) 
j=i 

Here, ipi nc is the incident wave, iPe(t, fj) is the ex- 
citing field for scatterer j, and T(fj) its scattering 
operator. Focusing on scatterer j = 1 as represen- 
tative, we have 

N 

tpE(f, n) = tp inc (f) + Y T{rj)il! E {r, fj). (6) 

If, as is the case here, the positions are not known, 
it becomes necessary to assume some type of sta- 
tistical description. The number density n for 
N discrete scatterers with centers at fi, f 2 ■ ■ ., is 
n(fi, f 2 , . . . , ?tf) = Np(fi,f 2 , ...,f N ), where p is 
the probability density. Conditional densities, de- 
fined by fixing the position of one scatterer, satisfy 2 - 
n(f 2 , . . . ,f N \fi) = (N-l)p(f 2 , . . . ,f N \fi). For our 
purposes we need only the conditional number den- 
sity of the j th scatterer if a scatterer is known to 
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be at f\\ n(fj |fi) = (N — l)p(fj \fi), j ^= 1, satis- 
fying J" dfj n(fj |ri) = iV — 1. The analog of © is 
then an equation for the configurationally averaged 
field, (ipE(r\n)), 

(ipE(r\n)) = ip inc (f) 

d fj n(fj\n )T{fj)(^E (r\fj , fi)). (7) 



As noted by£, "The fact that the exciting field with 
one scatterer fixed is given in terms of the field 
with two scatterers fixed is the basic difficulty in 
multiple scattering". We adopt perhaps the sim- 
plest solution to this quandary, the QCA (quasi- 
crystalline approximation), under which assump- 
tion ([7| reduces to 

(ip E {f\fi)) = i/j inc {r) 

dfj nir^Tir^Eir^)). (8) 



The scattering operator for every scatterer is as- 
sumed to have translational invariance with 



T(6)e lkx = ^i n T n H^\kr)e l 



^^—e^-VfiO), (9) 
and therefore the Fourier series for / is 

f{6) = Y. T " eine - ( 10 ) 

n 

The assumed form of the operator in © means 
that each scatterer has the same scattering behav- 
ior. We consider plane wave incidence, ipi nc {f) = 
Ae lkx . The exciting potential (ip E (r,rj)) satisfies 
the Helmholtz equation and is regular function at 
the point fj, it can therefore be expressed 



(tl) E (f\fj)) = A„{fj) J„(kpj 



Hps) 



(11) 



where pj = f — fj, 9(pj) = aig(pj). Substituting 
from into the configurational average equation 
(|8]) and using the addition theorem for cylinder 
functions^, yields 

Anin) = + V(-i)f r„ +p / nifjln) 
„ Js+ 



(12) 



on the half-space S + = {x > 0}. The effective 
wavenumber £ defines the coherent wave according 
to the assumed functional form for each A„ : 



A n (f j )=i n A n e i ^. 



Use of the addition theorem for cylindrical func- 
tions then implies 

A n e* Xl = Ae^ + £ T n+p A n+p L p ($, (14) 
p 

with 

L p(0 = rP / dfjii(fj\fi) 
Js+ 

x H^ikr^y^^e 1 ^!, (15) 

where fji = fj - fi, 9(fji) = arg(r}i). 

We assume a modified form of the hole correc- 
tion with hole radius b, 

,~,~s J not 1 + Hfji)) , r 3 i>b, 

n ( r i\ r i) = S n (I 6 ) 
[0, Tji < b, 

where h represents the deviation from the constant 
background value. It satisfies h(r) — > as r —> oo, 
with the stronger condition 

lim R- 2 / h(r)rdr = 0. (17) 
R^co J b 

Following, e.g£, it may be shown that the integral 
([TBI reduces to 



L p (£) = 2im 



e-k 2 



k 2 



(18) 



where 



N P (&) = Zbf p (Zb)HW(kb) - kbJ p (£b)HW (kb), 

(19a) 

Mp(&)= I J P (^r)H^(kr)h(r)k 2 rdr. (19b) 

Jb 

The simple "hole correction" is ft, = and hence 
M p = 0; see Sec. IV. D of Ref£ for further discus- 
sion of the case h ^ 0. 

B. Matrix formulation 

Substituting from (JTSJ) into (fbfj) and equating to 
zero the coefficients of e^ Xl and e lkxi yields two 
equations, known as the Lorentz-Lorenz law and 
the extinction theorem respectively, 

A " + tJ^2 E M n - P {$)T P A P = 0, (20a) 



.4 



(13) 



P— — CO 



k{z-k) z 

X3 ' p— — oo 



Y2 t p A p = °> ( 20b ) 
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where (|20a[) is satisfied by all n £ Z and 



Af p m = N p ^b) + (^-l)M p (^b). 



(21) 



Equation (|20a[) is a homogeneous system which de- 
fines the effective wavenumber and the associated 
infinite eigenvector, while Eq. (|20b[) defines the 
amplitude of the eigenvector in terms of the exci- 
tation amplitude A. The identity A/"_ p (£) = A/p(£) 
follows from the known properties of Bcsscl and 
Hankcl functions. Equation (|20a[) with h = => 
M p = is equivalent to the system of equations 
studied by Linton-Martin^ (their equation (71)). 
We focus on the solution of Eq. (|20a|) . 

We introduce non-dimensional parameters y and 
e which depend on b, along with some related vec- 
tors and matrices: 



y = 


m 2 - {kbf, 


(22a) 


e = 


—AiriQb 2 , 


(22b) 


b = 


T^a, 


(22c) 


u = 


T^e, 


(22d) 


Q(y) = 


(kb)~ 2 T 1/2 QT 1/2 , 


(22e) 



where the vectors a and e, the diagonal matrix T, 
and the symmetric matrix Q, are defined 



a= (..., A-!, A , A u ...y, 
e =(..., 1, 1, 1,...)*, 

_ f M m - n {jb) - l 



Q 



(23a) 
(23b) 
(23c) 

(23d) 



(C/fc) 2 - 1 ' 

Then Eq. (|20a[) can be expressed 

{y(l-eQ(y))-euu t }b = 0. (24) 

Setting the determinant of the matrix {..} of this 
infinite system of equations yields the desired dis- 
persion relation for £. We next reduce this to a 
simple and transparent form, and obtain an ex- 
pression for the amplitude vector a. 

C. Implicit solutions for £ and a 



Multiply (|24|) from the left by the inverse of (I - 
eQ) , yielding 



(l + wu')b = 0, 



(25) 



with infinite vector w = — ey _1 (I — eQ) _1 u. Not- 
ing that det(I + wu*) = 1 + w'u, we deduce that 
the solution for y can be expressed implicitly as 



y 



eQ(»)]" 



(26) 



The dependence upon b in Eq. (|26p may be re- 
moved by introducing alternative non-dimensional 
scalars y = y/(kb) 2 and e = e/(kb) 2 , that is, 



1, 



-Ainok 



(27) 



and writing Q = Q(y), in terms of which Eq. (|26l) 
becomes 



y 



e'KeT)- 1 -^®] X e. 



(28) 



This formula clearly splits the dependence upon 
the scattering matrix T, from that of multiple in- 
teractions, Q. The choice of the symmetric ma- 
trix Q, and hence Q, is motivated by the obser- 
vation that at leading order the elements Af n are 
equal: J\f n (kb) = 2 /(in) for all n. Despite the ap- 
parent pole at £ = k in Eq. (|23dp . the matrix 
Q(y) is a regular function of y at the origin since 
the limit and its derivatives exist as £ — > k. The 
solution of Eq. (|20a[) may be expressed in implicit 
form through cither of the identities (|26|) or (|28l) . 
We will find both useful in different circumstances. 
The latter is simpler for considering general prop- 
erties, while the former is useful for the particular 
limit as b — » 0. 

The identity 1 + w*u = also implies that the 
null vector of (f2"5)l is of the form b = aw (a ^ 0). 
The precise value of a follows from the extinction 
theorem (|20b[) . and the vector of amplitudes A n 
can then be expressed as 



a = 



2k 
k + l 



[(i-eQ(y)T)] \. 



(29) 



D. Asymptotic expansion 

We seek an asymptotic expansion of y = y e in 
powers of the small parameter e: 



eyi + e 2 y 2 + e 3 y 3 + 



The individual terms follow from (f2l)|) as 

1 r\ n ( 

Vn=^* — [e[l-eQ(y e )Y 



u. 



(30) 



(31) 



Using the expansion (1 — x) 1 = l + x + x 2 + ... 
in (|31|) and noting that the derivative is evaluated 
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at e = means that only a finite number of terms 
are necessary for a given n. Thus, 



1 H™ 

2/™ = 3 U *X~H £™Q™~ 1 (0) + e n ~ 1 Q n ~ 2 (eyi) 



n\ d e 
e 2 Q(«/i 



2 y n -2) 



(32) 



=0 



We only need to expand each term in this finite 
series to obtain its 0(e") contribution. This is 
simple for the first term e"Q"(0), and for the sec- 
ond it is e"- 1 Q"(eyi) e"yi^Q™~ 2 (0) where 

^g(O) = ^9(y)\ y - - Subsequent terms become 
more complicated but the procedure for finding 
their contribution is straightforward. Thus, 

y„ = u t (Q"- 1 (0) + y 1 A Q »-2 (0) 



+ 



vl d 2 



d 



2 Ay 1 ay 



dy 



(33) 



c=0 



The following expansion includes all terms up to 
fourth order in the small parameter e, 

eu*u + e 2 u*Q u 



y 



(34) 



+ e 3 [u i Q 2 ,u+(u*u)u t Q u] 
+ e 4 [u'Qgu + (u t Q u)u'Q[ ) u 

+ 2(u*u)u t Q Q u + -(u*u) 2 u*Qju] + . . . , 

where 

Q = Q(0), Q o = Q'(0), Q o ' = Q"(0). (35) 

The symmetry of Q* = Q has been used to simplify 
terms in the expansion (|34l) . 



E. Asymptotic expansion for finite kb 

Expanding the function Af p (£b) for small (£6 — 
kb) yields 

Qomn = D^_ n (kb), Q' 0mn = D^_ n (kb), (36) 
where 



VK 
T 



x 2 J p (x)H^'(x) + 2M p {kb) 



D£){x) = - l -D$\x) + l ^M v {kb) 



(37a) 



+ \ [p 2 - x 2 - i2nx i J p (x)H^{x)] . (37b) 



Higher derivatives may be found using the identity 

m 2 N; + {ib)N' p + {m 2 - P 2 )N p 

= -2^b) 2 J p ^b)H^(kb). (38) 

Using the above results, the expansion to 0(e 3 ) 
may be rewritten in terms of the areal density of 
scatterers, no, as 

£ 2 = k 2 + dino + d 2 n 2 + d 3 nl + . . . , (39) 

where 

di = -4»5^T B = -4i/(0) J (40a) 

n 

rf2 = -p E D ( °l n (kb)T m T n , (40b) 



641 ^ ^L(kb)D^ 
,4 ' Dni-n(kb)T m T n . 



rf 3 = fcT E D ^-n(kb)D^i p {kb)T m T n T p 

ni,n.p 



16di 



(40c) 



The dependence on fc6 and the scattering function 
is contained in the coefficients cq, d2, etc. 



F. Small fc& limit 



In this limit we derive the terms in the series 
£ 2 = k 2 + Smo + 5 2 n 2 + S 3 nl + 64^ + .... (41) 



It follows immediately from Eq. (|40|) that Si = 
d\, while 82 may be found by letting kb — > and 



using D ( p \x) = i£!+0(z) and D ( p >{x) = f - 

D p °\x)+0(x) as x — > 0. It is easier, however, 
to begin with the small kb expansion of Q. This 
allows us to deduce not only the terms in (|39|) to 
O(n-o) but the next one. Higher order terms can 
be found using the procedure described next. As 
kb — >• we have, using (|19a|) , 

•w-s (!)"'• < 42 » 

Hence Q becomes independent of &, 

Q ro „(y)^((l + y)"-l)/y, (43) 

from which it follows that 
1, 



Qomn = -j\m - n\, 

Q'omn = gl™ - n \(\ m - n \ - 2 ); 



(44) 



)L = ~ n|(|m ~ n| ~ 2)(|m ~ n| ~ 4) ' 
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etc. Now express the asymptotic expansion (f3"4"|) 
in terms of the nondimensional parameters y and 
e which do not depend upon b, 

y = etrT + e 2 e*TQ Te 
+ e 3 [e*T(Q T) 2 e + (tr T)e*TQ Te] 
+ e 4 [e*T(Q T) 3 e + (e*TQ Te)e t TQoTe 

+ 2(tr T)e*TQ TQ Te + - (tr T)V TQ 'Te] 

(45) 



O e 5 



The matrices Q , Q , etc. are defined in the same 
way as in (|35f for the matrix Q(y) of (|23d[) . The 
coefficients in Eq. (|4T|) then follow from eqs. (|44| 
and (SSI) as 



Si 



k 2 



Y^\m-n\ 



(46a) 
(46b) 



m,n,p 

) \m - n) T m T n 



k 4 



2fc 2 



(46c) 



32 



-g ^2 \m - n\\n - p\\p - q\T m T n T p T q 



«8<5i 



— ^ |m - n|(n - p) 2 T m T n T p 



7n,n,p 



E 



^1 i i3 2(52 / \2 

6fc 4 2fc 2 /c 2 ' 



T T 



(46d) 



G. The Linton-Martin formula generalized 

The coefficients in (HBJ) depend upon the far-field 
scattering function through its Fourier coefficients. 
We now show that the coefficients can be expressed 
in terms of the function / itself rather than its 
Fourier series. Thus, 

5 1 = -4i/(0), (47a) 
^ = ^[^^)± [mf (-0)i (47b) 



6* = 



m 
^k 4 



d(9cot- / d0cot-S( 



251 d2 [f(0)f(-e)] 9=o 5A 



k 4 d(9 2 
32 



2k 2 



(47c) 



^4 = 577 / d6»cot- / d6»cot- / d(9cot- 

7r 3 fc 6 y 2 J 2 J 2 





1=0 



61 



3wk 6 
25 2 d 2 
+ fc 4 d<9 2 

where 
S(M) 



£d9cot( e -)-^[f(e)f(-9)} 



\f((t\f( f)W 6lh 62 6163 



2 • 

(47d) 



1 d 2 



A 8989 

f{0) [f(S)f(S -9)- f(9)f(-9 - 9)] 
f(-9) [f(9)f(9 - 9) - f(-9)f(9 + 9)]}, 
1 8 3 



1(6,6,9) 

8 898989 

f(0)f(9 - 9) [f(-9)f(9 - 9) - f(9)f(-9 - 9)] 
+ f(0).f(-9 - 9) [f(9)f(9 - 9) - f(-§)f(9 + 9)] 
+ f(~e)f(6 + 9) [f(9)f(-9 - 9) - f (-§)/(§ - §)] 

+ f(~0)f(6 - 9) [f(-6)f(6 + 9) - f(9)f(9 - 9)] 



We now justify these expressions. 

The first identity for 8\ in (|47[) is obvious from 
the definition of Eq. (fl"U)) . Regarding S 2 , we note 
that the product of f(8) and f(—9) may be written 



f(6)f(-6)= £ E T n T s e 



i(n—s)9 



(48) 



n— — 00 S— — 00 



Interchanging the indices, implies that the double 
sum is 

00 00 

f(9)f(-9)= E T n T s cos(n-s)9. (49) 

n— — 00 s— — 00 

This identity is equivalent to Eq. (83) of£ but 
without the restriction T n = T_„ that was as- 
sumed there. The formula for S 2 follows from the 
relational (Eq. 3.612(7)) 



1 f 7 * 9 

— / d#cot - sin m9 = sgn(m). 
n Jo 2 



(50) 
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Turning to the coefficient £3 in (|46c|) . define 

s(«^ = -EEE( m - n )( n -p) 

m n p 
X T m T n T p e *( m_7l ) e e *( n_ P)^ 

' -<?)/(-§)]■ (51) 



<90<90 



The individual terms in the triple sum are of the 
form e 1 '" 6 , but an identical sum with the preferred 
dependence sin • • ■ 9 in each term can be obtained 
by an appropriate permutation of 2 2 = 4 expres- 
sions. The correct combination is 

^{S{6, 8) - S{0, -0) - S(-9, 9) + S(-9, -9)}, 

which reduces to 8(0, 9) as given. It may then be 
deduced, again using the identity (|50|) . that 



\m - n\\n - p\T m T n T p 



— TT / 



d0cot-S(0,0). (52) 



The second term in (|47cj) follows directly from (|49|) . 
Regarding the coefficient 64, the third term in the 
right member of (|47d[) follows in the same manner 
as the integral for 62, based on Eq. (|49|) . The 
second term in (|47d[) may be deduced using the 
relation 

a ^-EEEf^ 11 )^-?) 2 

m n p 



89 



x T m T n T p sin(m — n)0cos(n — p)9, (53) 



evaluated at 9 = in combination with the integral 
identity (f50]>. Note that 



90 k 



/"(o)^[/(0)/(-0)] 



/'(O)^[/(0)/'(-0) + /(-0)/'(0)] 
/(O)^[/'(0)/'(-0)] + [/(*)/("*)] 



Finally the first term in the right member of ()47dl) 
may be obtained using the same type of argument 
used for the identity (|52|) . The function analogous 
to S now has three arguments, 

^^>- ! EEEE( m - n )( fi -rf 

m n p q 
X(p- q)T m T n T p T q e i(m-n)8 e i(n-p)S e i(p-q)B 



8 3 



898989 



[f(0)f(9 -0)f (9 -8)f (-§)]. (54) 



and 2 3 = 8 permutations are required in order to 
arrive at the correct quadruple summation, yield- 
ing T(0,0,0 = ). 



III. SIMPLIFICATIONS BASED ON THE FINITE 
RANK OF Q 

The infinite matrix Q is in practice well approx- 
imated by a matrix of finite rank. This follows 
from the fact that the far-field scattering function 
f(9) is an entire function of the angular argument 9 
considered as a complex variable^. The scattering 
operator (matrix) T is therefore compact and has 
only a finite number of eigenvalues of finite size. 
At low frequency only the first few elements T m 
for m near zero are important (monopolc, dipole, 
etc.). Furthermore, as we will see in this section, 
Q is of rank 2 in the high frequency limit. There- 
fore, at any finite frequency the infinite system of 
equations is really not so in practice, and may be 
replaced by a finite system. We first develop the 
solution for finite rank n, and then apply it to two 
important cases; the low frequency Rayleigh limit 
(n = 3) and the high frequency limit (n = 2). For 
the remainder of the paper we take h, introduced 
in Eq. (fT6)h to be zero, so that M p = 0. 



A. Q is of rank n 

As noted above only a finite number of the ele- 
ments T n are significant at any given frequency. If 
only n are non-zero, then Q is of rank n. The ma- 
trix satisfies an homogeneous equation of degree 
n + 1, 



Q n+1 + ^Q n + 6 n -iQ r 



+ 6!Q = 0, (55) 



with 61 = (-l) n (detQ), b n = -(trQ), and 
therefore, 

(I - eQ)- 1 = I + lC Q + ^ 2 e 2 Q 2 + . . . + /3 n e n Q n , 
where 

1 + he n+1 - k 



1 ' Z^m=l u ™ fc 



(56) 



Thus, using the fundamental result (|26[) it fol- 
lows that y = (£ 2 — fc 2 )6 2 is given by 

n 

y = eu'u + J2 < ; i.yi'Q 1 "- (57) 
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Since Q = Q(y), (|57p is an implicit equation for y 
which could be solved by iteration, for instance. It 
may also be written as 



e = k 2 - 4m /(0) 



(58) 



1. Rayleigh limit, Q is of rank 3 

An important case is n = 3 which is useful at 
low frequency (Rayleigh limit) when the scattering 
matrix is well approximated by only three terms: 
To and T±\. For n = 3, 

Q 4 -IqQ 3 +IqQ 2 -I[qQ = 0, (59) 

with I Q = trQ, JL Q = J[(trQ) 2 -trQ 2 ], W Q = 
dct Q. In this case the implicit equation for y be- 
comes 

y = eu'u + e 2 [(l - el Q + e 2 l Q ) u*Qu 

+ (1 - elQ)eu'Q 2 u + e 2 u*Q 3 u] 

X [l-do + e'lg-^mg]" 1 . (60) 

More detailed results for the Rayleigh limit are pre- 
sented in ffV] for the particular case of circularly 
cylindrical scatterers. 



B. High frequency limit 

In the high frequency limit kb^b 3> 1, we have, 
from (|19a[) . 



-N p (&) tt^., A+(-lfB, (61) 



with 



A = (k+J) jtk-Qb B = ( fc ~J) e »(fe+C)b 



The high frequency limit for elastic waves was dis- 
cussed in£, and the same methods developed there 
could be used for the acoustic problem. It is in- 
structive to note that (foTj) and (|22e[) . combined 
with M p = 0, implies that 



where 



Q = auu 1 + /3vv', 

a=(A-l)/y, P = B/y, 

v = T 1 / 2 (..., -1, 1, -1„ 1, . . .)*, 



(62) 



1/2 

with vq = T . The matrix Q is therefore rank 
2, and the wavenumber £ can be found using the 
methods described above. 

Equation (|60|) reduces for rank 2 to 



t , 2 (l-etrQ) u t Qu + e u t Q 2 u 

?/ = euu + e — tt~, — ^ ■ 63 

y 1 - etrQ + e 2 dctQ v ; 

Noting that 

u*u = v'v = /(0), u*v = v*u = f(w), (64) 
and hence 

trQ = (a + /3)/(0), 
dctQ = a/3(/ 2 (0)-/ 2 W), 

the equation for £ becomes 

y = e/(0) + 6 2 x 

(a/ 2 (0) + /3/ 2 (^) - e/(0)a/3(/ 2 (0) - f 2 ^))} 
/{l - c /(0)(a + P) + e 2 a/3(/ 2 (0) - / 2 (tt)) }. 

This simplifies to 

(£ 2 -^) 2 +4ino(f 2 -fc 2 )(A + fl)/(0) 
-16n 2 AB(/ 2 (0)-/ 2 (7r)) = 0, (65) 

or dividing out the factor (£ 2 — k 2 ) (corresponding 
to the trivial solution £ 2 = k 2 ) implies that £ at 
high frequency satisfies 

e = k 2 - ^/(0) [(* + 0e"« 6 + i« - *)e* 6 ] e* fcfc 



+ ^(/ 2 (0)-/ 2 W)e^. 



(66) 



IV. GENERALIZED SELF CONSISTENT MODEL 



The Generalized Self Consistent Model (GSCM) 
developed by Yang and MaU2 was derived using 
a self consistent scheme applied to the Waterman 
and Truell's formula 2 -. Among the objectives of 
this section is to apply this scheme, which is very 
broad in scope, to the generalized Linton and Mar- 
tin formula. 

The basic idea follows Christensen and Lo^. In- 
stead of considering cylinders that are directly im- 
mersed in fluid, a "three phase cylinder" model 
is used which assumes that each cylinder is sur- 
rounded by a cylindrical ring of fluid, the whole 
being immersed in a outer region of equivalent fluid 
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of unlimited extent. Hence, form functions in this 
section correspond to three phase cylinders, and, 
as for cylinders, they may be expressed as a modal 
sum and calculated numerically^. 

Let a be the radius of circular cylinders and c 
(0 < c < 1) their concentration, the radius a c of 
cylindrical rings is related to a and c by 



n na 
ngna 2 



(67) 



Let p e f f and fc e / / be the effective properties of the 
equivalent outer fluid. The mass density /? e // is 
defined as the spatial average 

Peff = c Pcyl + (1 - c)Pfluid, (68) 

with p cy i and p fluid the cylinder and fluid mass 
densities. The wave number keff, which is un- 
known, is determined with the use of the self con- 
sistent scheme. 

Without loss of generality, consider the Linton 
and Martin formula at the second order in concen- 
tration. Let £lm be the Linton and Martin's ef- 
fective wave number in the outer equivalent fluid, 
then we have 



Clm = k lff ~ 4m /(fc e //,0) + 



8"o 
nk eff 



d6 



A d 



x ™t(-) — [f(k eff , 6) f(k eff , -6)}, (69) 



with 



eff) 



(70) 



and2£ no = c/ira 2 . The self consistent scheme con- 
sists in assuming that k e f / = £lm ■ From a physi- 
cal point of view this means that the outer equiva- 
lent fluid is a medium in which the waves propagate 
in exactly the same manner as the coherent waves. 
Because £lm = keff there is no scattering due to 
the three phase cylinders in the outer equivalent 
fluid, and the medium can be considered as ho- 
mogenized. It follows from (|69p that k e ff is given 
by the equation 



f(keff,0) 



2n 



d9 cot(-) 



X JqU '(kef f,6)f (k ef f, 



(71) 



It is worth mentioning that at low concentration 
of cylinders the second term in Eq. (|71|) can be 
neglected, so that Eq. (f7Tj) reduces to 



which corresponds to the equation for the CPA-i2. 
In other words, the CPA appears as the approxi- 
mation of Eq. (|TTj) to first order in concentration. 

Another way of presenting the self consistent 
scheme is to use an iterative procedure applied to 
Eq. (|69j). Starting with fco = k we carry out the 
homogenization by employing Eq. (|69l) to obtain k\ , 
and so on, so that 



k 2 



f(k ef f,0) = 0, 



(72) 



i -4m /(fc„,0) 
+ ^ J" d8 cot( -)±lf(k n ,O)f(k n ,-e)}. 

The iteration is repeated until there is convergence, 
k n+ i — > k e ff, the solution of (fTTj) . This proce- 
dure is of interest not only for computations but 
also for its physical interpretation. The effective 
wave number k\ corresponds to a coherent wave 
that accounts for the double interactions between 
cylinders according to the basic hypothesis of the 
Quasi Crystalline Approximation (QCA)2. In the 
same way, &2 involves double interactions in the 
medium characterized by k\ , the medium in which 
double interactions arc already taken into account. 
So, relative to the fluid medium characterized by 
ko, k% must deal with interactions of higher order, 
even if we do not know how to describe these ex- 
plicitly. In addition, the number of interactions 
between cylinders is linked to the order of conver- 
gence of the scheme. Accordingly, it can be sup- 
posed that the use of the self consistent scheme 
applied to an explicit effective wavenumber, what- 
ever it is, kisA, kwT or k^M (whatever the order in 
concentration), can improve the accuracy of the re- 
sults while the concentration of cylinders increases. 
This is what has been observed ir>i^ when search- 
ing for dynamic effective mechanical properties of 
composites at low frequencies. Furthermore, just 
as the Linton and Martin approach is more effi- 
cient than the one of Waterman and Truell, Eq. 
(jTTj) should be more efficient than the GSCM in- 
troduced by Yang and Mai, and also the CPA. Nu- 
merical calculations are beyond the scope of this 
paper, but the different effective wavenumbers are 
compared at low frequency in the next section. 



V. THE RAYLEIGH LIMIT 



In the Rayleigh or low- frequency limit, the size 
of the scatterers is assumed to be small compared 
to the incident wavelength. In this section scatter- 
ers are specifically considered to be circular fluid 



Multiple scattering by cylinders 9 



cylinders of radius a, characterized by the density 
P2 and the sound speed c 2 . They are immersed 
in a fluid, characterized by the density p\ and the 
sound speed c\. The Rayleigh limit then corre- 
sponds to k\a — > with k\ = u>/c±. It is then suf- 
ficient to take only the lowest order coefficients in 
the modal series^. More exactly, it can be shown 
that at low frequencies, only To and T±\ make a 
contribution, reducing the infinite matrix Q to a 
rank 3 matrix, see ijllll 

The goal of this section is to compare the 
Rayleigh limit for the different effective wave num- 
bers which correspond respectively to ISA, Wa- 
terman and Trucll, Linton and Martin, CPA, the 
GSCM developed by Yang and MaU^, noted here 
G-WT, and the GSCM introduced in JEV]from the 
Linton and Martin approach, noted here G-LM. 



A. Effective wave numbers 

We first calculate effective wave numbers using 
the Rayleigh limit of the far-field scattering func- 
tion associated with cylinders which are directly 
immersed in fluid: 



T 



±1 



iir 

T 

«7T / p 2 



P2C r 1)( " ia) ' 



pi 



92 + Pi 



(73a) 
(73b) 



These coefficients can be derived from those cor- 
responding to circular elastic cylinders immersed 
in fiuid3i, if the shear and longitudinal velocities 
are equal to and c 2 respectively (Scattering co- 
efficients are denoted by R in Reb^i with Ro = Tq 
and R±i = —T±i). 

The use of Eqs. (fT3")) leads to the following re- 
sults for the ISA, Waterman and Truell (WT) and 
Linton and Martin (LM). First 



fci 



= f 



Pic\ 



( k eff\2 



P2 ~ ipi 
Pi + P2 

(ISA), 



fci > ISA 

in terms of which the other two are 



(74) 



fci 1 



eff\2 



X < 



/ISA 
gig? _ 

P2cf 
2 

Pl£i _ 



2c' 



P2 ~ Pi 
P2 + Pi 



_2pi 



P2+P1 



(WT), 
(LM). 



(75) 



Note that the latter two effective wavenumbers are 
almost the same if the densities p\ and pi are close 
to one another in value. However, even at low 
frequency where cylinders look like "point scatter- 
ers" , the Waterman and Truell and the Linton and 
Martin approaches give different results as soon as 
p2 is not close to p\. 

B. Wave numbers from the self consistent scheme 

We now consider effective wave numbers ob- 
tained by the self consistent scheme, which uses 
scattering coefficients for the "three phase cylin- 
der" as described in SjlVI The coefficients are cal- 
culated as outlined ir>22, with the results 



Tn 



in 

T 



(l-c) 



Peff 
Pi 



, PeffC 2 l 
' P2C\ 



k 2 
fci 2 



T±i = —F(p eff )(k eff a c ) 2 , 



(kia c ) 2 , 
(76) 



where p e tf = (I — c)pi + cp2 from Eq. and 



(1 



F(Peff) = 
' c)(Pi - P2Peff) + (1 + C)p 1 (p 2 



Peff) 



(1 - c)(Pl + P2Peff) + (1 + C)P1{P2 + Peff) ' 

Note that Eqs. (J73J) follow from Eqs. (J7BJ formally 
if we first put k e ff = fci and p e // = Pi, in order 
to identify the outer equivalent fluid with the fluid 
itself, and then cancel the ring of fluid by setting 
a c = a, which implies c = I. 

We now consider three approaches based on self 
consistent schemes: the CPA, the generalized self 
consistent method based on Waterman and Truell 
(G-WT), and the same for Linton and Martin's 
approach (G-LM). We find 



eff< 



Peff 



ki 



[1 



Pi 

1 - ^(Peff 
1 + mPeff)} 
l+2F(p ef/ ) 



'P2C\ 



-2F2(p e// ) 



-l)cj 

(CPA), 
(G-WT), 

(G-LM). 



(77) 



In order to compare these effective wave numbers 
with those of eqs. (jT4"f and (jT5"j) in the previous 
subsection, we have to perform an asymptotic ex- 
pansion with regard to the concentration of scat- 
terers. It follows that 



F(peff) 



(Pi ~ P2? 

2pi 



Pi + P2 



Cf 

2p~i 
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and wc find at the second order in concentration 



fci 



x < 



pigi 



rpi£l 



2pi 



2c^ 



2 P2- Pi 



Pi 



P2+PI 

2pi 
P2+P1 

2pi - 
P2+P1 - 



(P1-P2) 2 I 
4pi(pi+P2) J 
(P1-P2) 2 1 
4pi(pi+p 2 ) J 



"Pi, 

(CPA), 

(G-WT), (78) 
(G-LM). 



As expected, all three methods give the same re- 
sult at the first order in concentration. At the sec- 
ond order, the CPA and G-WT introduce the same 
additional term as compared to Linton and Mar- 
tin's approach, but with the opposite sign in each. 
It is significant that the self consistent scheme 
applied to Linton and Martin's formula does not 
modify the result, at least at this order. The same 
cannot be said of the CPA and the Waterman and 
Trucll methods. Thus, we may conclude that Lin- 
ton and Martin's approach and the G-LM can be 
considered "self consistent" methods. Finally, we 
note that the additional terms in the CPA and in 
the G-WT in Eq. (fTSf are very small if the den- 
sities pi and P2 are close in value. Hence, all the 
methods are equivalent at low frequency when the 
densities p\ and P2 are equal. Of course, these re- 
sults say nothing about what happens at higher 
frequency and at higher concentration. 



VI. CONCLUSIONS 

Implications of the quasi-crystalline approxima- 
tion (QCA) on the effective wave number beyond 
the dilute concentration limit have been described. 
Equation ([25)1 is the starting point for all further 
developments, and as such represents the funda- 
mental result of the paper. It splits the implicit 
form of the effective wavenumber into two distinct 
parts, one defined by the single scatter T- matrix, 
T, and the other by the spatial arrangement of 
the scatterers, Q. In this paper we have used to 
the hole correction, for which Q is given by (|23d[) . 
More generally, this matrix has elements 



i(k-$)xi 



1 



2fc(£ - k) £ 2 - K 



(79) 



where £ m _ n is defined by Eq. (|15[) for arbitrary 
pair correlation function. The QCA is exact for 
a regular array of scatterers, in which case Q mn 
can be reduced to a known lattice sum. Equation 
(|28[) therefore provides a formula for determining 
the dispersion curves of a regular array. This and 
other implications will be examined elsewhere. 
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